Nonparametric assessment of mangrove ecosystem in the context of coastal resilience in Ghana

Abstract Cloud cover effects make it difficult to evaluate the mangrove ecosystem in tropical locations using solely optical satellite data. Therefore, it is essential to conduct a more precise evaluation using data from several sources and appropriate models in order to manage the mangrove ecosystem as effectively as feasible. In this study, the status of the mangrove ecosystem and its potential contribution to coastal resilience were evaluated using the Google Earth Engine (GEE) and the InVEST model. The GEE was used to map changes in mangrove and other land cover types for the years 2009 and 2019 by integrating both optical and radar data. The quantity allocation disagreement index (QADI) was used to assess the classification accuracy. Mangrove height and aboveground biomass density were estimated using GEE by extracting their values from radar image clipped with a digital elevation model and mangrove vector file. A universal allometric equation that relates canopy height to aboveground biomass was applied. The InVEST model was used to calculate a hazard index of every 250 m of the shoreline with and without mangrove ecosystem. Our result showed that about 16.9% and 21% of mangrove and other vegetation cover were lost between 2009 and 2019. However, water body and bare land/built‐up areas increased by 7% and 45%, respectively. The overall accuracy of 2009 and 2019 classifications was 99.6% (QADI = 0.00794) and 99.1% (QADI = 0.00529), respectively. Mangrove height and aboveground biomass generally decreased from 12.7 to 6.3 m and from 105 to 88 Mg/ha on average. The vulnerability index showed that 23%, 51% and 26% of the coastal segment in the presence of mangrove fall under very low/low, moderate and high risks, respectively. Whereas in the absence of mangrove, 8%, 38%, 39% and 15% fall under low, moderate, high and very high‐risk zones, respectively. This study will among other things help the stakeholders in coastal management and marine spatial planning to identify the need to focus on conservation practices.


| INTRODUC TI ON
For natural capital accounting, structural factors such as height of mangrove canopy and biomass density are important considerations.
They also provide soft defence for flood protection and play significant roles in controlling the global carbon balance to combat climate change (Houghton, 2005). Globally, deforestation of mangrove forest has been noted as a major source of greenhouse gas emissions (IPCC, 2007) and as a danger to the ecological services that mangroves supply (Aja et al., 2022). The evaluation of plans intending to lessen the effects of climate change using natural solutions requires accurate quantification and monitoring of temporal and spatial variations in mangrove cover and biomass (Omar et al., 2014). Fatoyinbo and Simard (2012) estimated that mangroves cover about 7600 ha along the coast of Ghana and support coastal populations who are at risk of disasters by providing important ecosystem services such as natural coastline protection, coastal erosion reduction, improvement of water quality and provision of alternative livelihoods (Aja et al., 2022;Mondal et al., 2018).
The ecosystem services provided by mangroves in Ghana are essential to the well-being of both local communities and the broader environment. Protecting and restoring these important ecosystems is therefore crucial for sustainable development in the region.
Mangrove assessment in the context of coastal resilience involves evaluating the role of mangrove ecosystems in promoting the ability of coastal areas to withstand and recover from natural disasters and other environmental stresses (Sebesvari et al., 2019;Whelchel et al., 2018). This assessment typically involves a range of factors, including the extent and diversity of mangrove ecosystems, the degree of protection they provide against storm surges and erosion and their ability to buffer against other climate-related impacts (Estrella & Saalismaa, 2013;Sudmeier-Rieux et al., 2021). To assess the role of mangroves in promoting coastal resilience, researchers have used a variety of methods, including remote sensing, field surveys and modelling.
It has been underlined how crucial remote sensing is in this context as one of the key sources of spatial data (Lagomasino et al., 2016). Given the complexity of the tropical mangrove forest ecosystem, remote sensing methods are not without their limits when assessing mangrove ecosystem (Son et al., 2015). The tropical environment's complex mangrove forest ecosystem, cloud cover and saturation at particular biomass levels make it difficult to use optical satellite technology only for mangrove monitoring and assessment (Aja et al., 2022). The closed canopy structure of tropical mangrove forests and the fact that spectral signature is insensitive to the spatial changes of aboveground biomass (AGB) greater than 150 Mg/ha may be the reason for the lack of effectiveness of optical remote sensing in mangrove structural assessment (Omar et al., 2014).
Aboveground biomass is estimated differently through a backscatter mechanism using active remote sensing techniques such as synthetic aperture radar (SAR). The ability of SAR technology to penetrate through closed tree canopy has contributed to the recent advances in mangrove assessment and monitoring (Omar et al., 2014).
As SAR is one of the practical ways to acquire remote sensing data without the effect of cloud cover and independent of weather or light conditions, interest in radar remote sensing for monitoring tropical mangrove forests is developing (Thomas et al., 2017). LiDAR is another method with the ability to precisely capture the structure of mangroves without weather interference (Tianyu et al., 2020), but the lack of time series data with global coverage limits its application on a continental to global scale (Hu et al., 2016).
Many researchers are now incorporating radar data from various sensors for a more accurate quantification of mangrove extent and other biophysical properties due to the recent availability of active satellite data with increasing spatiotemporal coverage (De Santiago et al., 2013;Pinki et al., 2019;Thomas et al., 2018). The majority of satellite estimations of AGB were generated from the shuttle radar topographic mission digital elevation model (SRTM DEM) that was produced in the year 2000 (Fatoyinbo and Simard, 2013;Giri et al., 2011). The digital elevation model can be combined with synthetic aperture radar images, such as ALOS PALSAR data, to further improve the satellite estimations of aboveground biomass density (Omar et al., 2014).
Conventional mapping methods have been employed by many researchers to provide information on mangrove structure, but these methods are constrained by the availability of images, the requisite computational power and the required technical know-how (Gorelick et al., 2017;Yancho et al., 2020). In particular for mangrove mapping and monitoring, new techniques and technology such as cloud computing are ushering in a new era (Wang et al., 2019;Wulder et al., 2018;Yancho et al., 2020). Platforms for geospatial data processing on the cloud, such as Google Earth Engine (GEE), give users unparalleled access to a vast library of ready-to-use geospatial data (Gorelick et al., 2017;Yancho et al., 2020).
Google Earth Engine eliminates the requirement for local standalone computers to download and process data by storing data locally and performing client-versus-server operations (Amani et al., 2020;Yancho et al., 2020). This resolves issues with hardware and the requirement for technical know-how necessary for remote sensing. These advancements enable the rapid development and application of different mapping techniques across various spatial scales, opening up advanced remote sensing applications to more users (Gorelick et al., 2017;Hansen et al., 2013). Additionally, GEE open-source nature allows for flexibility and customization while facilitating methodological repeatability (Vos et al., 2019). Currently, a number of studies have used GEE to map mangrove AGB and shown its efficiency, offering encouraging advancements over earlier techniques. Carreiras et al. (2013) used L-band synthetic aperture radar coupled with field data to estimate AGB and reduced uncertainties. Fatoyinbo and Simard (2013) used a combination of mangrove map derived from Landsat, lidar canopy height estimates and elevation data from SRTM to derive the structure and biomass of mangroves for the African continent.
Assessing coastal vulnerabilities and the potential for mangrove ecosystem to contribute towards the protection of coastal communities is essential for effective long-term planning, sustainability and the management of coastal resilience (Al Ruheili & Boluwade, 2023). The Integrated Valuation of Ecosystem Service and Trade-offs (InVEST) Coastal Vulnerability model has been applied elsewhere to assess coastal vulnerability and to explore the role of natural habitats such as mangroves to mitigate coastal hazards (Arkema et al., 2013;Aneseyee et al., 2022;Ai et al., 2022;Al Ruheili & Boluwade, 2023). Mangrove ecosystems must be accurately mapped and inventoried in order to determine how they are changing and to what pressures they are subjected to, before their potential role as a natural defence can be assessed. Such information serves as the foundation for in-depth social-ecological assessments that can guide efforts to lessen anthropogenic pressures on mangrove and encourage sustainability and conservation.
Few researchers have combined optical and radar data for improved mapping and assessment of mangroves (Ghorbanian et al., 2021;Yancho et al., 2020); however, such research as well as information on the role of natural habitats in protecting coastal areas is limited in Ghana. In this study, we used GEE to combine optical and radar data set in a spatial framework to monitor and evaluate mangrove ecosystem. We employed InVEST model to simulate the degree to which mangrove ecosystem contributes to coastal resilience. The aim of this article was to analyse and evaluate mangrove ecosystem with a view to determine the impact of human-induced pressure on the temporal and spatial variation of mangrove extent and AGB density as well as to demonstrate the importance of the mangrove ecosystem as a robust nature-based solution to coastal disasters. The research questions that this study intends to address are: (1) What is the impact of human-induced pressure on the temporal and spatial variation of mangrove extent and aboveground biomass (AGB) density? and (2) How does this affect the ability of mangrove ecosystems to serve as nature-based solutions to coastal disasters? 2 | MATERIAL S AND ME THODS

| Description of study site
This study was conducted at a complex wetland fringing part of Anlo community, Shama district, Western region, Ghana ( Figure 1).
The 52 km 2 delineated area is bounded within latitude 5°1′30″ N and 5°3′5″ N, longitudes 1°34′30″ W and 1°37′30″ W. The wetland is situated at the lower part of the Pra watershed with an area of 23,000 km 2 (Opoku-Ankomah, 2000) and occupies roughly 20% of the landmass of Ghana (Owusu et al., 2015). In terms of hydrology, this area is within the floodplains of Pra River, which has a direct link with the ocean.
The study area witnesses yearly high temperature ranging from 21.74°C to 31.6°C (Bessah et al., 2018) with a monthly average relative humidity above 70% (Kwabena et al., 2013) and a long-term average rainfall of 1446 mm (Bessah et al., 2019), making the climatic condition a tropical monsoon (Am) (Kottek et al., 2006). The dominant soil types in the upper part of the basin are Acrisols and Lixisols, which are characterized by fairly high possibility for run-off (Ross et al., 2018), exposing the low parts of the watershed to water disasters during periods of heavy rainfall (Bessah et al., 2020).
The elevation is very low (0-42 m), the coastline is irregular with a sandy beach next to the wetland and the ocean is marked The community heavily depends on the mangrove ecosystem for F I G U R E 1 Study location map. a variety of activities, which intensifies the dynamics of land use change (Aja et al., 2022).
htm. These data are a dual-polarized L-band SAR with a 25-m resolution, 10 0 × 10 0 in longitude and latitude. The image collection ID is N06W002_09_sl_HH, N06W002_09_sl_HV. The Google Earth Engine database, which is accessible at https://devel opers.
google.com/earth -engin e/guide s/senti nel1, was used to retrieve the Sentinel-1 data for 2019. These data are a C-Band synthetic aperture radar with a descending pass, dual-polarized with a 25 -m resolution in interferometric wide swath mode. ee.ImageCollection ('COPERNICUS/S1_GRD') is the Image Collection ID. The database of the GEE was used to retrieve the 2009 Landsat 7 surface reflectance tier1 data. Four visible and near-infrared bands, two short wave infrared bands and one thermal infrared band are all present in this data. ee.ImageCollection ('LANDSAT/LE07/C01/T1_SR') is the image collection ID: More information about this data can be found at https://www.usgs.gov/lands at-missi ons/lands at-surfa cerefle ctance.
Landsat 8 surface reflectance tier1 data (2019) were retrieve from GEE database. These data contain five visible and near-infrared bands, two short wave infrared bands and two thermal infrared bands. The image ID is ee.ImageCollection ('LANDSAT/LC08/C01/ T1_SR'). Additional information about these data can be found at https://www.usgs.gov/lands at-missi ons/lands at-surfa ce-refle ctance. The SRTM DEM data were retrieved from GEE database. in 2010 was needed to delineate mangroves and was accessed at https://data.unep-wcmc.org/datas ets/45 (Giri et al., 2011). Global mangrove watch is a shapefile of global mangrove areas as at 2010.
We used bathymetry from the general bathymetric chart of the global ocean (GEBCO) available at https://www.gebco.net/data, wavewatch III and continental shelf contour to assess coastal vulnerability. The GEBCO bathymetry has a spatial resolution of 15 arc sec. Information about wavewatch III and continental shelf contour can be found at Coastal Vulnerability Model-InVEST documentation (stora ge.googl eapis.com).

| Methods
For this work, mangrove evaluation and monitoring were made using Google Earth Engine (GEE), a cloud-based technology. To study the spatial changes in mangrove cover over time, LULC maps were made by jointly classifying SAR and optical data for two periods, namely 2009 and 2019. Mangrove aboveground biomass (AGB) as at year 2000 was estimated using basal area-weighted height estimates from the SRTM DEM (Simard et al., 2019). The AGB estimates from satellite was compared with field measurements to examine the temporal and spatial variation over time.

| Mangrove extent mapping
The method involves retrieving the following data set: Sentinel-1 data, ALOS PALSAR-1 and 2 data, Landsat 8 and Landsat 7 surface reflectance tier1 data, and the mangrove distribution vector from GEE database as shown in Figure 2. These data were loaded into the Google Earth Engine code editor and then these data were preprocessed using cloud masking for the optical images (Giri et al., 2011;Hansen et al., 2013) and SAR speckle filtering to reduce speckle noise (Ayman et al., 2017). The images were enhanced by creating a function that masks out cloud shadows from optical images and another function was created that filters the speckle noise from SAR images.
The ALOS PALSAR-1 and 2 data were converted from log 10 to dB (decibels) using the formula below (JAXA, 2006).
The NDVI (normalized difference vegetation index), which served as an information layer for the classifier, was used to transform the optical data into a composite image (Shi et al., 2016;Thomas et al., 2018).

| Random Forest (RF) algorithm and land use land cover change analysis
An ensembled tree-based machine learning algorithm called random forest (RF) that employs decision trees to perform image classification was used. In this supervised classification, the algorithm builds a decision tree for each sample according to the predictors, the trees cast votes for each pixel to identify the type of land cover and the most supported value is given to each pixel. In order to 'train' the RF classifier, training samples must be created (Pelletier et al., 2016;Shelestov et al., 2017).
The methodology employed in this work is consistent with that suggested by Ghorbanian et al. (2021) and Yancho et al. (2020), and it entails gathering backscatter samples that represent each relevant landcover class and using the samples to train the classifier. With 100 trees and five predictors chosen at random for each split, the RF algorithm was run (Ghorbanian et al., 2021). Field campaigns that took place between December 2020 and April 2021 served as the basis for the training and validation data used in this work. To help with the creation of 'training' samples and result validation, high-resolution The output was transferred to google drive to be processed further.

| Classification accuracy and result validation
The overall accuracy and the quantity allocation disagreement index (QADI) were calculated using the following formula: where, N, total number of sampled pixels; A, allocation disagreement of classes; Q, quantity disagreement of classes.
The quality allocation disagreement index (QADI) classification accuracy scale is as shown in Table 1. The quality allocation disagreement index is an efficient classification accuracy assessment for a variety of machine learning algorithms including random forest (RF). It has advantage over the traditional kappa coefficient because it solves some of the problematic paradoxes associated with Kappa (Feizizadeh et al., 2022).
Overall, 2131 training sample points were produced. The sample points were separated at random into 'training' and 'validation' groups of which 80 per cent of the points were used to 'train' the model, while 20 per cent were used for validation (Mahdianpari et al., 2020). As a result of using the same pixels for classifier training and validation, any systematic mistake was eliminated by this process (Geiß et al., 2017;Pimple et al., 2018;Shi et al., 2016). Utilizing 20% of the training sample, stratified random samples were created for an independent accuracy evaluation, which follows the procedure outlined by Barenblitt and Fatoyinbo (2020). High-resolution satellite imagery that is accessible in ArcGIS was used to verify each point.

| Estimation of mangrove stand height
Mangrove canopy height is a key element in calculating aboveground biomass and carbon sequestration rates since it has a strong correlation with carbon turnover through litterfall Overall accuracy = Number of correctly classified pixels (sum of diagonal) Number of total sampled pixel × 100 (3) production (Rovai et al., 2018;Saenger & Snedaker, 1993). In this study, mangrove canopy height estimate was made using SRTM digital elevation measurement (Farr et al., 2007) generated in February 2000. Mangrove height estimation was made following an approach that had had been effectively applied at regional scales (Fatoyinbo and Simard, 2013;Simard et al., 2006Simard et al., , 2008. The SRTM DEM values represent a height that lies between the elevation of the ground and the canopy's highest point (aka Lorey's height) (Lagomasino et al., 2016;Simard et al., 2006Simard et al., , 2008. The maximum height is estimated to be 1.697 × (times) this value based on empirical data from field measurements (Simard et al., 2019).
We used the global mangrove distribution vector to isolate mangrove areas and mask nonmangrove regions in the SRTM elevation data set (Giri et al., 2011).
To estimate mangrove height, the SAR image (ALOS PALSAR-1) and the SRTM elevation model were clipped to the global mangrove distribution vector file ( Figure 2

| Mangrove aboveground biomass estimation
Aboveground biomass (AGB) estimates can be produced using allometric equations once the height of the mangrove canopy has been determined (Comley & McGuinness, 2005). The present study estimates AGB based on the correlation between canopy height and mangrove biomass (Lucas et al., 2014;Simard et al., 2019). We where H ba , basal area-weighted height (~1.08*SRTM).

| Field inventory
Mangrove biomass is measured either directly through destructive harvesting of specific trees or indirectly through measurements of tree diameters and inference using allometric relationships (Chave et al., 2005). The destructive approach requires cutting down the trees, but the nondestructive approach uses allometric equations (Gibbs et al., 2007). Allometric equations, however, frequently depend on the locality and tree species ( Table 2 because none of the indicated wood densities are particular to Ghana.

| Coastal vulnerability assessment using InVEST model
In order to compare the relative exposure of the coastlines with and without mangrove ecosystems, we estimated the coastal exposure

| Accuracy assessment
To offer thorough statistical analysis for each classification, confusion matrix and an independent accuracy assessment were used.
The user and producer accuracy were quite high according to the confusion matrix for the two time periods (Tables 3 and 4). However, according to the independent accuracy assessment, somewhat higher accuracy was discovered for 2009, which made use of ALOS PALSAR-2 and Landsat-7 (Figures 7 and 8).

| Estimation of canopy height
Tree height, density and basal area are often taken into account when evaluating the structure of a forest; these factors are also required to calculate the amount of aboveground biomass. We created a mangrove canopy height map for the year 2000 using the global mangrove distribution vector file, SRTM DEM, and SAR data. This was then utilized to create estimations of AGB for the same time period. Our analysis of mangrove canopy height distribution shows that in year 2000, mangrove canopy height ranges from about 6.3 to 12.7 m (Figure 9).
This estimate of mangrove height significantly corresponds to previously reported values in different estuaries in Ghana (Aheto et al., 2011 andNortey et al., 2016). The local-scale geophysical determinants of environmental gradients, such as nutrient availability and soil or water salinity, are highlighted by the spatial variability in canopy height (Thom, 1982). The spatial distribution of aboveground biomass patterns is determined by the relative influence of local environmental gradients, as well as the spatial distribution of mangroves species in a given coastal region (Rovai et al., 2018).
The result of AGB estimate for year 2000 showed that aboveground biomass ranged from 0 to 105 Mg/ha ( Figure 10).   (Table 5). Aboveground biomass ranged from 0 to 87.52 Mg/ha. The AGB estimates generally decreased with distance from the shoreline.

| Coastal vulnerability assessment
We calculate the index of vulnerability for every 250 m along the coast in our study area. This index was used to identify the role of mangrove ecosystem as a nature-based solution for coastal defence and resilience. The result is presented in Figures 11 and 12, the areas in purple colour denote high-risk areas, whereas areas in yellow/white colour denote low-risk areas. Figure 11 showed that in the presence of Mangrove ecosystem, 23% of the coastal segment fall under very low-/low-risk zone, 51% of the coastal segment fall under moderate risk, while 26% fall under high-risk zone. Figure 12 showed that in the absence of Mangrove ecosystem, 8% of the coastal segment fall under low-risk zone, 38% fall under moderate-risk zone, 39% fall under high-risk zone, while 15% fall under very high-risk zone.
F I G U R E 6 Map of land cover extent for 2019.

| DISCUSS ION
The importance of combining optical and radar data for long-term coastal ecosystem mapping was highlighted by Mwita et al. (2012), Wijedasa et al. (2012) and Aja et al. (2022). Combining radar and optical imagery allows for effective landcover mapping, and recent developments in artificial intelligence, machine learning algorithms and the advent of platforms for high-performance cloud computing, such as Google Earth Engine (GEE) have made this possible (Erika et al., 2020;Midekisa et al., 2017;Wijedasa et al., 2012). This gave us the opportunity to thoroughly evaluate the ecosystem of interest.
The To offer thorough statistical analysis for each classification, we made use of confusion matrices and an independent accuracy evaluation. The confusion matrix for the 2009 imagery classification revealed that, of the 647 pixels that were identified as mangroves, all were classified correctly, whereas the confusion matrix for the 2019 image revealed that, of the 647 pixels that were identified as mangroves, 642 pixels were classified correctly as mangroves (Tables 3 and 4). This could be as a result of the effect of using L-band F I G U R E 7 Independent accuracy assessment for 2009 land cover classification.
F I G U R E 8 Independent accuracy assessment for 2019 land cover classification.
in 2009 and C-band in 2019 classification. Similarly, Attarchi and Gloaguen (2014) as well as Aja et al. (2022) noted that L-band synthetic aperture radar performs better than C-band when combined with optical data. The strength of the classification method is underscored by the independent accuracy assessment (Figures 7 and 8 risk areas represented 23% of the coastal segment, the moderate (rank = 3) risk areas represented 51%, whereas the high-risk areas (rank = 4) represented 26% (Figure 11). In the absence of mangrove ecosystem as shown in Figure 12, very low-/low-risk areas (rank = 1 & 2) represented 8%, moderate-risk areas (rank = 3) represented 38%, high-risk areas (rank = 4) represented 39%, while very high-risk areas (rank = 5) represented 15%. As shown in Figures 11 and 12, coastal vulnerability will increase if mangrove forest is removed. This result agrees with the findings of Arkema et al. (2013) and Al Ruheili and Boluwade (2023).

| CON CLUS ION
The rate of tropical and subtropical mangrove forest degradation is increasing, and this is increasing the chance that coastal disasters will harm more people and property. Consequently, there is a critical need to map and keep track of this ecosystem and to identify their multiple benefits and values in the context of coastal resilience. In this study, we mapped and inventoried a complex mangrove wetland based on satellite data from multiple sources, secondary data sets, state of the art modelling tools and appraised the impacts of anthropogenic activities on their temporal and spatial changes. By so doing, we identified that the mangrove ecosystem is being converted to other land uses, which is also leading to a decline in the AGB density and these patterns of change can help guide conservation efforts and enable sustainable usage of our priceless mangrove ecosystem.
Knowing where mangroves are most likely to lessen vulnerability to erosion and floods from storms and future sea levels, as well as safeguard vulnerable people and property, is necessary to prioritize mangrove habitats for conservation or restoration in service of natural hazard reduction.
We developed an index of hazard that takes into account the protective function of mangrove ecosystems at every 250-m scale in order to pinpoint the coastline segments where mangroves have the greatest potential to safeguard coastal communities against hazard. We evaluated the vulnerability of the coastline with and without mangrove ecosystem for the entire coastline of 593 points along the coastal segment. We found that in the presence of mangrove, about 74% of the coastal segment fall under low-/moderate-risk zone, which decreased to 46% without mangrove.
For prioritizing mangroves for coastal defense, the vulnerability index from this paper is most helpful. Our findings highlighted that continuous degradation of existing mangrove ecosystems will lead to greater harm to people and property. We therefore recommend that coastal defense planning in Ghana must incorporate mangrove ecosystems restoration alongside physical structures.
The accuracy of our result is within the limited bounds of the data spatial resolution. A major limitation of this study is the lack of up-todate images (e.g. 2020 images) and high-resolution images (e.g. 10-m resolution) in GEE database for the study region as at the time of this analysis. The coastal vulnerability suit in the InVEST model serves as an effective tool for rapid assessment of habitat scenarios in the context of coastal resilience to inform decision-makers. However, the model has some limitations because of the use of a surge proxy that could lead to a generalization of storm dynamics and may result in the overestimation or underestimation of risks for some coastal areas. Future study will extend the methodology to the entire coastline of Ghana, incorporating the impacts of sea level rise.

ACK N OWLED G M ENTS
We are grateful for the free access to the datasets provided by Google earth engine developers and the Natural capital project. We would like to thank the Africa centre of excellence in coastal resilience for being the enabler. We thank the two anonymous expert reviewers for their constructive comments which helped to improve the quality of this work.

CO N FLI C T O F I NTER E S T S TATEM ENT
We have no conflicting interest to declare.

DATA AVA I L A B I L I T Y S TAT E M E N T
The links to the data sets analysed during this study are provided in